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In the framework of the fct-factorization approach, the production and polarization of T mesons at the 
Fermilab Tevatron is considered, and a comparision of the calculated double differential distributions and 
spin alignment parameter a with the DO experimental data is shown. We argue that measuring the double 
differential cross section and the polarization of upsilonium states can serve as a crucial test discriminating 
two competing theoretical approaches to the parton dynamics in QCD. 



Nowadays, the production of heavy quarkonium 
states at high energies is under intense theoretical and 
experimental study [H [2] . The production mechanism 
involves the physics of both short and long distances, 
and so, appeals to both perturbative and nonpertur- 
bative methods of QCD. This feature gives rise to two 
competing theoretical approaches known in the litera- 
ture as the color-singlet [3l|4] and color-octet [5] models. 
According to the color-singlet approach, the formation 
of a colorless final state takes place already at the level of 
,the hard partonic subprocess (which includes the emis- 
sion of hard gluons when necessary). In the color-octet 
model, also known as nonrelativistic QCD (NRQCD), 
the formation of a meson starts from a color-octet QQ 
pair and proceeds via the emission of soft nonpertur- 
,bative gluons. The former model has a well defined 
applicability range and has already demonstrated its 
predictive power in describing the J/ip production at 
HERA, both in the coUinear [B] and the /ct -factorization 
[7] approaches. As it was shown in the analysis of 
recent ZEUS [5j data, there is no need in the color- 
octet contribution, neither in the coUinear nor in the 
fct-factorization approach. The numerical estimates of 
the color octet contributions extracted from the analysis 
of Tevatron data are at odds with the HERA data, es- 
pecially as far as the inelasticity parameter z = / 
is concerned H]. In the /cj-factorization approach, the 
values of the color-octet contributions obtained as fits 
of the Tevatron data appear to be substantially smaller 



than the ones in the coUinear scheme, or even can be 
neglected at all [H [H [H US] . 

Recently, the results of new theoretical caclula- 
tions of the ncxt-to-leading (NLO) and next-to-next- 
to-leading (NNLO) order corrections to colour singlet 
(CS) quarkonium production have been obtained in the 
framework of standard pQCD [14] . In the region of mod- 
erate pt {pt > 10 GeV), these corrections enhance the 
color singlet production rate by one order of magnitude 
and even larger. These new results are in much better 
agreement with the fct-factorization predictions than it 
was seen for leading order coUinear calculations. 

In the present note we follow the guideline of our pre- 
vious publication |15| and show a more detailed analysis 
of the production and polarization of T mesons at the 
Tevatron conditions using the fcj-factorization approach. 

2. THEORETICAL FRAMEWORK 

In the fci-factorization approach, the cross section of 
a physical process is calculated as a convolution of the 
off-shell partonic cross section a and unintegrated par- 
ton distribustions JFg(a;, fc|,, /i^), which depend on both 
the longitudinal momentum fraction x and transverse 
momentum fc^: 



X (T(^(^(xi, X2 , k-^rp ^ k^rp ^ ■■■) dx'2 dk^z-p dk<2'p. (f) 

In accordance with the fct-factorization prescriptions 
O [III [iHl [19] , the off-shell gluon spin density matrix 
is taken in the form 
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In all other respects, our calculations follow the stan- 
dard Feynman rules. 

In order to estimate the degree of theoretical uncer- 
tainty connected with the choice of unintegrated gluon 
density, we use two different parametrizations, which are 
known to show the largest difference with each other, 
namely, the ones proposed in Refs. [16l[T9] and [20| . 

In the first case [16], the unintegrated gluon den- 
sity is derived from the ordinary (collinear) density 
G{x, /i^) by differentiating it with respect to /x^ and set- 
ting jj? = k^. Here we use the leading order Gliick- 
Reya-Vogt (LO GRV) set [31] as the input colinear 
density. In the following, this will be referred to as 
dGRV parametrisation. The other unintegrated gluon 
density [20] is obtained as a solution of leading order 
Balitsky-Fadin-Kuraev-Lipatov (BFKL) equation [19] 
in the double-logarithm approximation. Technically, it 
is calculated as a convolution of the ordinary gluon den- 
sity with some universal weight factor. In the following, 
this will be referred to as JB parametrisation. 

The production of T(liS') mesons in pp collisions can 
proceed via either direct gluon-gluon fusion or the pro- 
duction of P-wave states Xb followed by their radiative 
decays Xfc^f +7- The direct mechanism corresponds to 
the partonic subprocess g + g ^ + 5 which includes 
the emission of an additional hard gluon in the final 
state. The production of P-wave mesons is given by 
g + g ^ Xb: and there is no emission of any additional 
gluons. As we have already argued in our previous pub- 
lication [15], we see no need in taking the color-octet 
contributions into consideration. 

The polarization state of a vector meson is character- 
ized by the spin alignment parameter a which is defined 
as a function of any kinematic variable as 

a{P) = {da/dV - 3daL/dP)/{da/dP + daL/dV), (3) 

where a is the reaction cross section and ci, is the part of 
cross section corresponding to mesons with longitudinal 
polarization (zero helicity state). The limiting values 
a = 1 and a = —1 refer to the totally transverse and 
totally longitudinal polarizations. We will be interested 
in the behavior of a as a function of the T transverse 
momentum: V = \pt\- The experimental definition of 
a is based on measuring the angular distributions of the 
decay leptons 



dr(T^/i+^")/dcos( 
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where 9 is the polar angle of the final state muon mea- 
sured in the decaying meson rest frame. 

The definition of helicity and, consequently, the def- 
inition of a is frame-dependent. There are four com- 
monly used different definitions of the helicity frame: 



these arc the recoil, the target, the CoUins-Soper, and 
the Gottfried- Jackson systems. In our analysis, we will 
basically use the recoil system (which, at the Tevatron 
conditions, is the same as the laboratory or proton- 
proton center-of-mass system), unless a different choice 
is explicitly declared. 

When considering the polarization properties of 
T(IS') mesons originating from radiative decays of P- 
wave states, we rely upon the dominance of elec- 
tricdipole El transitions"^' . The corresponding invariant 
amplitudes can be written as [22 



iA{xi 
»-4(x2 



T7) oc 6^''"'3fc,,4'^^)em4''\ 



(5) 
(6) 



with p and k being the momenta of the decaying meson 
and the emitted photon; e'^^', e'^^^', e^^', and e^"*"' the 
respective polarization vectors; and e^"'"/' the antisym- 
metric Levita-Civita tensor. This leads to the following 
relations between the production cross sections for dif- 
ferent hehcity states (see Eq. (14) in JT?): 



crT(/i=0) = B^Xl^^l) 

+ S(X2^T7) 
ar{\h\=i) = -B(Xi^T7) 
+ S(X2^T7) 
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The dominance of electric dipole transitions (at least 
for the charmonium family) is supported by the recent 
experimental data collected by the E835 Collaboration 
[53] at the Fermilab. 

All the other essential parameters were taken as in 
our previous paper: the 6-quark mass mi, = mx/2 = 
4.75 GeV; the T meson wave function |5't(0)P = 0.4 
GeV"^ (known from the Icptonic decay width r/+/- [21]); 
the wave function of P-wave states |\l/^(0)p = 0.12 
GeV^ (taken from the potential model [25]); the radia- 
tive decay branchings Br{xb,j^^l) = 0.06, 0.35, 0.22 
for (J = 0, 1, 2) [24 ; the renormalization and factoriza- 
tion scale = ji^p = ji^ = m\ -I- p^. 



•^'in our previous paper 1151 . two somewhat difTcrent models 
were used for this process. 
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3. NUMERICAL RESULTS 

The results of our calculations are presented in Figs. 
1-4. Fig. 1 displays the px dependence of the differ- 
ential cross section and spin alignement parameter a 
for four different intervals of rapidity. Complementary 
to Fig. 1, Fig. 2 exhibits the rapidity dependence of 
the cross section and parameter alpha for three differ- 
ent intervals of px- Everywhere, we separately show 
the contribution from the direct production mechanism 
taken solely (thin curves) and after having the x decays 
added (thick curves). When possible, we compare our 
theoretical predictions with experimental measurements 
[26]-[28]. 

First of all, we notice the importance of the feed- 
down from Xb deccays, without which the experimen- 
tal data can hardly be understood. The calculations 
seem to underestimate the cross section data by ap- 
proximately a factor of 2. This can be considered as 
a room for higher order corrections and contributions 
from other possible subprocesses, such as the associated 
production oi T + b + b states. The latter was shown to 
be comparable in size with the ordinary production at 
high pt [29] . Any way, the disagreement by a factor of 2 
must not be taken too seriously,as it lies within the un- 
certainty connected with the choice of factorization and 
renormalization scales The JB gluon density leads 
to significantly better agreement with the data than the 
dGRV density. 

While the direct and indirect production mecha- 
nisms lead to more or less similar px and y spectra, 
the behavior of the polarization is very much different. 
This is seen in the right parts of Figs. 1 and 2, and is 
vividly shown in Fig. 3. 

Our results for the direct mechanism are also ap- 
plicable to the production of T(3S') states, with the 
only exception that the overall dimuon rate Br^^a{T) 
is lower by an approximate factor of 4 because of smaller 
value of the wave function (|4'(0)3sp : |\E'(0)isp oc 
ri+i-{3S) : r,+,-(lS') = 0.44 : 1.34 [24]) and smaller 
branching fraction (2.18% versus 2.48% [24]). In this 
case, the absense of the feed-down from Xb decays would 
make the experimental sample cleaner and clearer for 
theoretical analysis^-'. 

We also have to draw attention to the fact that the 
behavior of the spin alignement parameter a is frame 
depenent, as is demonstrated in Fig. 4. In particular, 
the sharp dip of a at y = is only seen in the recoil 



■''As it has been argued in |30| . the proper choice should be 
rather fi^ = (M^ +p^)/2 than +Px- 

^'Unless there exists an unobserved Xbi-iP) state, still below 
the open BB threshold, which is yet not excluded. 



system, but not in cither of the other three helicity sys- 
tems. This property has to be not forgotten in order 
that the comparison between the theoretical and exper- 
imental results be fully adequate. 

4. CONCLUSIONS 

We have considered the production of T mesons in 
high energy collisions in the /ct -factorization approach 
and compared the predictions on the differential cross 
sections and spin alignment parameter a with new DO 
and CDF data. Wc find a more or less reasonable agree- 
ment in all cases. 

We have argued that measuring the double differ- 
ential cross sections and, especially, the polarization of 
quarkonium states in extended px and rapidity inter- 
vals can provide interesting and important information 
on their production mechanisms. 

The purest probe is provided by the polarization of 
T(3S') mesons. In that case, the polarization is the 
strongest and the predictions are free from uncertain- 
ties coming from radiative Xb decays. 
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Differential cross section and spin alignement parame- 
ter a as functions of the T(IS') transverse momentum 
Pt, integrated over four different rapidity intervals. The 
panels from top to bottom: \y\ < 0.6; 0.6 < \y\ < 1.2; 
1.2 < \y\ < 1.8; 1.8 < \y\. Dashed histograms, dCRV 
gluon density; dash-dotted histograms, JB gluon den- 
sity. Thin lines, the direct contribution only; thick lines, 
with the feed-down from Xb states added. Experimental 
points; • DO [21]; o CDF [23; * DO (preliminary) [55] • 



Differential cross section and spin alignement parame- 
ter a as functions of the T(15') rapidity y, integrated 
over three different intervals of pr ■ The panels from top 
to bottom: pt < 3 GeV; 3 < pt < 8 GeV; 8 < pr 
GeV. Notation of the curves is as in Fig. 1. Recoil 
system is assumed everywhere. 



Fraction of longitudinally polarised T(IS') mesons 
ci(7(helicity=0)/dcr(all helicities) as function of the 
transverse momentum pr and rapidity y. Upper panel, 
direct subprocess; lower panel, Xb decays solely. 



Rapidity dependence of the parameter a as seen in the 
different helicity frames (sole Xb contribution with JB 
gluon densities). Dash-dotted histograms, recoil sys- 
tem; dashed histograms, target system (equivalent to 
Gottfried- Jackson system); dotted histograms, Collins- 
Soper system. 
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